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In a single-component superfluid under rotation a broken symmetry in the order parameter space 
results in a broken translational symmetry in real space: a vortex lattice. If translational symmetry 
is restored, the phase of the order parameter disorders and thus the broken symmetry in the order 
parameter space is also restored. We show that for Bose- Condensate mixtures in optical lattices 
with negative dissipationless drag, a new situation arises. This state is a modulated vortex liquid 
which breaks translational symmetry in the direction transverse to the rotation vector. 



An important property of a superfluid is its specific 
rotational response. Namely it comes into rotation by 
means of the formation of a vortex lattice. Under the in- 
fluence of other factors such as temperature, multiplicity 
of superfluid components, inhomogeneities etc., different 
"aggregate" states of vortex matter may form, such as 
vortex liquids, glasses, etc Q. The variety of states is 
even richer in multicomponent systems Q. The transi- 
tions between the various "aggregate" states of vortex 
matter are related to various ordering processess of par- 
ticles in condensed matter systems. For example, the 
process of thermal vortex lattice melting can be mapped 
onto an insulator-to-superfluid transition of bosons. In 
this mapping, a vortex line is viewed as a world line of 
a boson with the z-axis mapped onto a "time" -axis and 
the vortex liquid, which is also entangled, represents the 
delocalized/superfluid state of the dual bosons 3] . The 
central result of the present work is that we find evidence 
in large-scale Monte Carlo (MC) computations that in 
a two-component Bose-Einstein condensate (BEC), the 
vortex lines support a state possessing properties of a vor- 
tex liquid simultaneously with properties of vortex lattice 
i.e. breakdown of translational symmetry. 

Recent progress in creating and observing various mix- 
tures of multi-component BEC has produced much inter- 
est in these systems. When BEC components are not 
spatially separated, the generic type of interaction be- 
tween them is the current-current interaction (Andreev- 
Bashkin effect) describing non-dissipative drag be- 
tween the two superfluid components. Such a system in 
units where h = 1 is described by the free energy density 
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where mi, ra 2 , 0\ and 2 are the masses and the phases of 
the condensates while ni, n 2 control the phase stiffnesses 
of the two components. Further, the drag coeficient rid 



controls the density of one component dragged by the su- 
perfluid velocity of the other component. Since here we 
are interested in the physics of rotating system we include 
the field which accounts for rotation with angular ve- 
locity ft = V x = 27rfz, where mif is the number 
of rotation-induced vortices per unit area of component 
i. We will use / = 1/64 throughout. In the following, 
we denote vortices with 2ttI^^ windings in #(i,2) with a 
pair of integers (A6 1 = 2tt/i, A<9 2 = 2tt/ 2 ) = (h,h). The 
last term in (pQ) is the current-current interaction [4j which 
may be caused by different reasons, such as intercompo- 
nent van der Waals interaction 0] or originating with the 
underlying optical lattice @|. It was first considered in 
the contexts of the physics of 3 He - 4 He mixtures and co- 
existing neutronic and protonic condensates in neutron 
stars. Various aspects of the rotational response of this 
system has been studied so far only for positive values 
of drag in context of 3 He — A He mixtures [4] and BEC 
mixtures 0- However, it has been recently shown that 
in optical lattices there arises an intriguing possibility 
to produce a BEC mixture with a negative inter-species 
drag n d 



In this paper we address the physics of a rotating sys- 
tem with a negative intercomponent drag and find it be- 
ing very rich. This is manifested in the situations where 
the usual notions from disordered versus ordered vortex 
states do not directly apply. Let us first briefly recapit- 
ulate the phase diagram of this system in the absence of 
rotation. Its main feature is that for significantly large 
drag \nd\ > n c , the easiest topological defects to excite 
thermally are (1,1) vortex loops. Proliferation of these 
composite defects leads to a state with order only in the 
phase difference, a so-called super-counter-fiuid [1, 0, @]. 
In order to estimate the phase stiffness which is left in the 
system after the (1,1) vortex loops proliferate, one has to 
extract from the Eq. (1) the term which depends only 
on the gradients of the phase difference, which stiffness 
is thus unaffected by the proliferation of (1,1) loops. The 
corresponding separation of variables in the presence of 
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After proliferation of (1,1) vortices it is the first term 
in ([2]) which accounts for the only phase stiffness remain- 
ing in the system, and we can renormalize the coefficient 
of the second term to zero and discard it. The complex- 
ity of the situation arising under rotation is that along 
with the (1,1) vortex loops excitations there are rotation- 
induced vortex lines. Vortex loops and lines affect each 
other's orderings and proliferation. Thus, we may ask (i) 
what are the ordering patterns of rotation-induced vortex 
lines in the model (pQ), (ii) how do rotation-induced vortex 
lines contribute to renormalization of stiffness and thus 
to symmetry breakdown patterns, and (hi) can ordering 
of the rotation-induced vortices signal the presence of a 
negative drag effect. 

To address these questions, we have performed large- 
scale MC computations using discretization of Eq. (pQ) un- 
der rotation, in the Villain approximation [6]. Through- 
out, we use a temperature scale such that the temper- 
ature T at which two decoupled superfluids with equal 
masses and phase stiffnesses transition to normal fluids 
is T = 3.3. A negative intercomponent drag will tend to 
increase the temperature at which this transition occurs. 

Consider first the simplest limit when mi = m 2 = 1 
and n i = n 2 = n. Eq. ([2]) then simplifies to 
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When nd = 0, the condensates are decoupled and a rotat- 
ing system forms two hexagonal lattices of the types (1,0) 
and (0,1). For rid < 0, an attractive interaction between 
rotation- induced vortices results. Thus, in the ground 
state the system forms a triangular lattice of (1,1) vor- 
tices. Such a configuration minimizes the gradients in the 
first term in Eq. (|3|). In the simplest limit mi = m 2 = 1 
and 774 = Ti2 = n we have found regimes where the vor- 
tex lattice melts while vortices nonetheless retain their 
composite character. The system retains a symmetry 
in the phase difference thereby representing a "rotation- 
induced" super-counter-fluid state. Introducing a mass 



and density disparity mi ^ m 2 and ri\ ^ n 2 gives an en- 
tirely different ordering and symmetry breakdown. This 
is the main focus of this paper. 

We study the spatial symmetry breakdown pattern 
and the effect of thermal fluctuations, by computing real 
space averages of vortex densities. They are produced 
by integrating the z-directed vortex segments along the 
z-axis, Vi(r±) = L~ x v£(r±, z), with a subsequent 
averaging over typically 10 4 different configurations at a 
given temperature. v?(r±,z) is the vorticity of compo- 
nent i in the z-direction at r = (x, y, z) and r± = (x, y) 
is the position in the x/y-plane. Thus, for an elementary 
vortex on the numerical grid directed along z-axis, the 
quantity v£(r±,z) is nonzero and positive in lattice pla- 
quette which corresponds to the center of the vortex. It 
is nonzero and negative for an antivortex, whence Ui(r±) 
gives the average a^-position density of the rotation- 
induced vortices. 

Let us consider the case n 2 /ni = 4, n^/ni = —5.0 
and m 2 /mi = 2. Now, since the vortex density is 
proportional to flrrii there are twice as many vor- 
tices in component 2 as in component 1 for these param- 
eters. The system exhibits a striking vortex ordering. 
Component 1 forms a triangular lattice, while compo- 
nent 2 with twice as many vortices, forms a honeycomb 
lattice. Every second vortex in the honeycomb lattice is 
co-centered with a vortex of the other component. This 
can be also viewed as an ordered equal mixture of (1,1) 
and (0,1) vortices. We find that the structure with a 
honeycomb plus hexagonal vortex lattice persists for a 
significant range of temperatures. Fig. [T^i shows a real- 
space average and a 3d snapshot of a typical configu- 
ration of this spatial symmetry breakdown pattern at 
T = 6.9. This ordering has broken down at T = 9.5, 
where we observe a partial meltdown manifested in the 
disappearance of every second vortex position peak in the 
real-space averages. However, every other vorticity peak 
corresponding to a hexagonal sublattice co-centered with 
vortex lattice of component 1 survives. The reduction of 
the number of vorticity peaks in component 2, the corre- 
sponding change in the structure factor, along with a 3d 
snapshot of a typical vortex configuration, is shown on 
Fig. Ub- The structure factor of component i is defined 
as S«(fc_L) = \l/(L x L y f) Z r± Vi(r ± )e- [r ^ | 2 . 

Let us finally turn to the case where m 2 /mi = 2, but 
Ti2/ni = 16. Now, with n^/ni = —2.5, we find that at 
low temperatures, the system instead forms two square 
lattices. In the ground state the square lattice in compo- 
nent 2, which has twice as many vortices as component 1, 
is rotated 45 degrees with respect to the lattice of com- 
ponent 1, so that every second vortex is co-centered with 
a vortex of component 1 lattice, see Fig. [2^. Again this 
can be viewed as an equal mixture of (1,1) and (0,1) vor- 
tices. Note that, in contrast to the case of two-component 
vortex matter with only repulsive interactions @, OH, 
here the vortex lattices are not interlaced. Here, the 




FIG. 1: (color online) Vortex orderings for the parameters 
712/ ni — 4, rid/ni — —5.0 and 1712/ mi — 2. Panel (a) shows 
the low temperature phase, (T = 6.90), where vorticity aver- 
ages (shown in gold for component 2 and green for component 
1) produce a hexagonal vortex lattice in one component and 
a honeycomb lattice in the other component. Panel (b) shows 
the situation occurring at a higher temperature (T = 9.52), 
where a hexagonal sublattice of the honeycomb lattice melts, 
the number of vorticity peaks in real-space averages dimin- 
ishes by a factor two, and the remaining hexagonal lattice is 
co-centered with a vortex lattice of component 1. The the 
insets in the top panels in part (a) and (b) show the structure 
factor evolution. The lower panels show a typical 3d vortex 
configuration in a 24 x 24 x 24 segment of the system. 



FIG. 2: (color online) Vortex orderings for the parameters 
712/711 — 16, Ud/ni — —2.5 and 7712/7711 — 2. The fc-space 
inset is the full fc-space structure factor for component 1, in 
both panel (a) and (b). The left hand side in panel (a) and (b) 
is a realspace average of component 2 (yellow) , while the right 
hand side is component 1 (green). Column (a) corresponds 
to the ordering in low-temperature (T = 11.0). Coulumn (b) 
shows the situation taking place at higher temperature (T = 
13.3) where the green vortices (component 1) are in the phase 
corresponding to a dual supernuid bosonic density wave. It 
is far from obvious from the 3d picture in the panel (b) that 
component 1 breaks translational symmetry. However, when 
averaged over many configurations, we see that component 1 
indeed breaks translational symmetry. 



appearance of square symmetry is caused by attractive 
interactions between vortices of different types. 

As the temperature is increased from T = 11.0 to 
T = 13.3, cf. Fig. [21 the evolution of the system 
is particularly remarkable: we observe a discontinuous 
phase transition, where in the real-space averages the 
number of vortex position peaks in component 1 doubles. 
This should be compared with the previous case of lower 
disparity of stiffnesses. There, in contrast, the system 
undergoes a transition to a state where vortex position 
peaks of the other component was reduced by a factor 
of two. Furthermore, both lattices change symmetry by 
collapsing onto a hexagonal co-centered configuration, as 
seen in the right panel of Fig. [21 A 3d snapshot of a 
part of the system, shown in the lower panel in Fig. [21 
reveals that the process is accompanied by a rapid in- 
crease of vortex loops in component 1. Furthermore, the 
figure [3] shows the central feature of this state. Namely, 
the helicity modulus, equivalently the supernuid density 
computed according to the procedures in Ref. 0, for 
component 1 disappears essentially simultaneously with 
the structure factor for the square lattice in component 
1. However, at the same time there emerges a nonzero 
triangular structure factor in component 1. It extends 
for a significant range of temperature where the helicity 
modulus of component 1 is zero. Therefore, the above 



observations are not related to a standard vortex-loop 
proliferation transition in the 3dXY universality class. 
If this behavior were associated with a standard vortex- 
loop proliferation transition of the vortices in component 
1, the superfluid stiffness (helicity modulus) of this com- 
ponent would vanish simultaneously with the structure 
factor of the corresponding vortex lattice. 

Thus, we have a quite remarkable situation. On the 
one hand, zero helicity modulus in z-direction indicates 
that vortices are entangled with each other like in a vor- 
tex liquid, a state which has a dual counterpart in super- 
fluid bosons [3j. On the other hand, the vortex system 
nonetheless features a structure function which is char- 
acteristic of a vortex lattice, namely it has distinct peaks 
at reciprocal lattice vectors. Thus, the dual counterpart 
of the vortex system we found, is a bosonic superfluid 
density wave. 

In terms of vortex matter this corresponds to the fol- 
lowing situation. Vortices in component 1 in this state 
are largely co-centered with the vortex lattice of compo- 
nent 2, but at these temperatures constantly and freely 
switch from being co-centered with one to being co- 
centered with the another vortex at different points along 
the z-axis. In order for the number of vortex position 
peaks of component 1 to be double the number that is 
generated by the rotation, a large number of (1,0) vor- 
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tex loops must be induced. Then the part of each (1,0) 
loop that is parallel to the (0,1) vortices have a tendency 
to be co-centered with a (0,1) vortex, breaking transla- 
tional symmetry for this segment, while the remaining 
part of the (1,0) loop which is not parallel to the (0,1) 
vortices has a random position and does not break trans- 
lational symmetry. Only at further elevated tempera- 
tures, a crossover takes place where vortex loops prolifer- 
ate, and the vortices loose line tension and the structure 
factor vanishes, see Fig. El 
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FIG. 3: (color online) Plot of the evolution of the struc- 
ture factor for component 1 (red and green) when with in- 
creasing temperature the system undergoes a sequence of 
the phase transitions from the state with U(l) x U(l) bro- 
ken symmetry to the state with U(l) broken symmetry and 
finally to the fully symmetric state (denoted as "Sym"). 
The (green) squares represent the structure function of the 
square lattice in component 1, Si\\ at the Bragg-peaks at 
k ± = (0.7854,0.00) « (7r/4,0). The (red) crosses repre- 
sent the structure function of the triangular lattice in com- 
ponent 1, S£], at the Bragg-peak at k ± = (1.1781, 0.0982) & 
(37r/8, 7r/64). Furthermore, the (blue) circles represent the 
helicity modulus in the ^-direction for component 1. It van- 
ishes at the same temperature as the square lattice ordering 
ceases. However at the same temperature there appears a 
nonzero structure factor for a triangular lattice S^l • This is 
accociated with the vortex state dual to bosonic superfluid 
density wave. The parameters are n^l 'n\ — 16, n^/ni = —2.5 
and rri2/mi — 2. The system size is L x L x L, with L — 64. 
Periodic boundary conditions are used in all directions, 10 5 
sweeps are used for thermalization and 10 6 sweeps for collect- 
ing average values with sampling every 100 th sweep. 

In conclusion, we have considered two-component su- 
perfluids with a negative dissipationless drag. In the 
model Eq. (pQ) , the underlying optical lattice plays only a 
microscopic role by providing a negative intercomponent 
drag through the mechanisms discussed in Thus , 

Eq. (pQ) describes the system at temperatures T larger 
than the vortex pinning energy of optical lattice E p and 



thus there is no lattice pinning effect [ll| . The vortex or- 
dering pattern in this system is strongly affected by the 
negative dissipationless drag, resulting in the formation 
of square and honeycomb lattices. Observation of these 
different ordering symmetries in experiments would be 
the hallmark of intercomponent drag. At finite tempera- 
ture there are phase transitions between states with dif- 



ferent lattice symmetries. The main conclusion of our pa- 
per is that, apart from different patterns of spacial sym- 
metry breakdown, the standard notions of vortex order- 
ing single-component vortex matter do not directly apply 
in the case of two-component vortex matter with a neg- 
ative drag. Namely, we have identified a state of vortex 
matter which is dual to a bosonic superfluid density wave, 
where one of the components breaks translational sym- 
metry even though there is no symmetry broken in the 
order parameter space. In this regime, a standard exper- 
imental technique of a density snapshot with a significant 
averaging along the z-axis would indicate a vortex lattice 
even though this is not a superfluid state. Since this state 
is phase-disordered, it can be discriminated from super- 
fluid vortex lattice via interference experiments. In an 
experimental situation, these effects will be naturally af- 
fected by density inhomogeneities present in traps. How- 
ever, studies [12| of the effect of the presence of a trap on 
three dimensional vortex matter, suggest that the above 
states can be realizable in an extended area near the cen- 
ter of the trap. 
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